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ABSTRACT 

Context. It is widely accepted that the onset of the explosive carbon burning in the core of a carbon-oxygen white dwarf (CO WD) 
triggers the ignition of a type la supernova (SN la). The features of the ignition are among the few free parameters of the SN la 
explosion theory. 

Aims. We explore the role for the ignition process of two different issues: firstly, the ignition is studied in WD models coming from 
different accretion histories. Secondly, we estimate how a different reaction rate for C-buming can affect the ignition. 
Methods. Two-dimensional hydrodynamical simulations of temperature perturbations in the WD core ("bubbles") are performed 
with the FLASH code. In order to evaluate the impact of the C-burning reaction rate on the WD model, the evolution code 
FLASH_THE.T0RT01SE from Lesaffre et al. (2006) is used. 

Results. In different WD models a key role is played by the different gravitational acceleration in the progenitor's core. As a con- 
sequence, the ignition is disfavored at a large distance from the WD center in models with a larger central density, resulting from 
the evolution of initially more massive progenitors. Changes in the C reaction rate at T < 5 x 10** K slightly influence the ignition 
density in the WD core, while the ignition temperature is almost unaffected. Recent measurements of new resonances in the C-burning 
reaction rate (Spillane et al. 2007) do not affect the core conditions of the WD significantly. 

Conclusions. This simple analysis, performed on the features of the temperature perturbations in the WD core, should be extended 
in the framework of the state-of-the-art numerical tools for studying the turbulent convection and ignition in the WD core. Future 
measurements of the C-buming reactions cross section at low energy, though certainly useful, are not expected to affect dramatically 
our current understanding of the ignition process. 

Key words, supernovae: general - hydrodynamics - methods: numerical - nuclear reactions, nucleosynthesis, abundances - white 
dwarfs 



1. Introduction 

The sequence of events that leads to the ignition of ther- 
monuclear flames in the core of mass-accreting carbon-oxygen 
. white dwarfs (CO WD) marks the start of the explosion 
' of a type la supernova (SN la). It is known from theory 
; dffillebrandt & Niemevedl2000l) that the ignition process spans 
■ over a broad range of typical length and time scales, correspond- 
' ing to the transition from hydrostatic to explosive C-burning. As 
, a consequence, a proper modeling of this evolutionary phase is 
very challenging both from the theoretical and computational 
point of view. In the framework of the single degenerate explo- 
sion scenario, a detailed knowledge of the features of the WD 
at ignition is crucial for determining the initial conditions of the 
explosion models, i.e. for shaping the initial flame morphology 
in the numerical simulations. 

The dependence of the explosion outcome and features on 
the initial spatial distribution of the flame kernels has been 
often debated, given the underlying relative freedom in con- 
staining the unknowns of the multi-spot ignition scenario 
JOarc fa-Senz & Bravo 2005; Ropke et al. 2006b, 2007b). In 
general, it turns out that a larger number of ignition points re- 
sults in a more vigorous deflagration phase. Interestingly, in re- 



cent models with deflagration to detonation transition (DDT), 
ignition closer to the center leads to a weaker detonation phase 
and a smaller production of ^^Ni, becaus e of the stronger WD 
pre-expansion in the deflagration stage (IWooslev et al.1 120081 : 
iKasenet al.ll2009h . 

Several different approaches have been used to address the 
problem of SN la ignition. Most of the work on this topic has 
been devoted to full-star si mulations of the CO WD, both us- 
ing analytical models ( Wooslev et al.l l2004t I W unsch & Wooslevl 
200 4i) and numerical simulation s (Hoflich & Stein 2002; 
Kuhlen et al.'l2006t JPiro & Bildstenll2008»: ,Piro & Chana, 20081: 
Zingale et al. 200^. A parallel branch of studies focuses 
on the role of the convective Urea process for the ignition 
(.Lesaffre et al...2005a.b, and references therein). 

Recently iLesaffre et al.l (l2006h studied the link between the 
evolution of the progenitor in its binary system and the igni- 
tion conditions of SNe la. They found a broad range of densi- 
ties and temperatures for WD models at ignition, linked to the 
initial WD mass and the binary accretion history. The effect of 
the progenitor diversity on the explosion phase still needs to be 
explored. In particular, a fundamental issue has to be fully ad- 
dressed, namely the relation between the properties of the igni- 
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tion process and the homogeneity of SNe la as a class of objects. 
Does ignition contribute to the homogeneity of SNe la or, on the 
contrary (cf. Ropke et al. 2006a), is it more intrinsically linked 
to the residual diversity of Branch-normal SNe la? 

A complementary contribution to the study of this prob- 
lem is provided by the small-scale (1 km or less) modelling 
appro ach. These m odels are motivated by analytical calcula- 
tions (iGarcfa-Sen z & Woosley 1995; Wunsch & Woosley 2004) 
which predict that temperature perturbations ("bubbles"), pro- 
duced in the WD's core by the turbulent convective flow, finally 
trigger the thermonuclear runaway. The physics of the therma l 
fluctuati on has bee n addressed by IGarcia-Senz & Bravol (l2005h 
and Zin gale & Dursi (2007) . which focus on the transition from 
hydrostatic to explosi ve burning and the pr opagation of the first 
thermonuclear flames. 'lapi chino et alj (1200 6) studied the bubble 
physics with a series of small-scale 2D simulations. According 
to their results, the main feature of the bubble evolution is the 
fragmentation caused by its own buoyant rise motion, which can 
be in a sense considered as a cooling effect, counteracting the 
heating from the hydrostatic nuclear burning. Among nuclear 
burning and hydrodynamical dispersion, the process with the 
smaller timescale determines the outcome of the bubble, namely 
if it will trigger a thermonuclear runaway or will be dispersed 
and cooled down by thermal conduction. 

The goal of the small-scale studies is the description of the 
SNe la ignition scenario in terms of number of flame kernels, 
size of the ignition region in the WD core, (a)symmetry of the 
ignition and stochastic behaviors of the process. There is general 
agreement that the WD convective flow prior to the thermonu- 
clear runaway plays a crucial role on the ignition. It has also been 
proposed that the flame kernels could tend to cluster spatially, 
thus evolving in a sort of ignition plume, rather than separate 
flame seeds (e.g.. lLivne et al. I I2OOI . This idea is closely linked 
to the flow in th e WD convective core, which is Ukely to develop 
a dipole f eature (Wooslev et al. '2004': IWunsch & Wooslevl2004l: 
iKuhlen et al...2006: Zingale et al .,20091) ! 

The study of llapichino et al.l (l2006h explored only the igni- 
tion parameters that are inherent to the bubble physics: the initial 
bubble diameter, the initial bubble temperature, the initial dis- 
tance of the bubble from the WD center. Though instructive for 
understanding the small-scale physics of the temperature per- 
turbations, that work could only very generically constrain the 
large-scale ignition features, through a number of additional the- 
oretical assumptions. In the present work, we widen the scope to 
further ignition parameters, which can be related to the ignition 
condition in a more straightfor ward way. In particu lar, the bub- 
ble simulations developed in lap ichino et al.l (l2006l) will be used 
to probe the i gnition conditions in the WD models presented in 
iLesaflre et al.i (2006). We will show how the large range of igni- 
tion densities in the aforementioned WD models can be linked 
to a significant diversity for the size of the ignition region. 

In the second part of the present work, these results will be 
applied to a largely unaccounted issue of the ignition process, 
namely the assessment of the role of experimenta l uncertainties 
in the C-burning reaction rates. In a recent studv. ISpillane et aP 
(1200 7*) present new reaction rates for the reactions '^C('^C, 
a)^''Ne a nd ^^C('^C, p)^^'Na , which are larger than previously 
assumed dCaughlan & Fowleiill988h . Those authors suggest that 
the new measurements could affect the ignition and even the ex- 
plosion theory of SNe la. We show in this work that the impact of 
these measurements on the WD evolution to ignition is actually 
margin al, and by using the WD evolution code by Le saffre et a l. 
(I2OO6I) we make estimates on the role of a hypothetical larger (or 
lower) rate. 



Table 1. Main features of the four WD models used in the sim- 
ulations. 



Model 


T, [10** K] 


[lO-'gcm-^^] 


g [10^ cm s -] 


rb [km] 


1 


7.84 


1.82 


5.18 


163 


3 


7.63 


2.08 


5.90 


154 


4 


7.15 


3.12 


8.78 


135 


6 


6.94 


4.17 


11.7 


123 



The structure of the paper is the following: in Sect. |2] the 
simulations will be introduced, with special emphasis on the de- 
tails related to the WD background models (Sect. 12.2b and the 
C-burning reaction network (Sect. l231 l. The result are presented 
in Sect. [3] and the conclusions are drawn in Sect.|4] 



2. Numerical tools 

2.1. Bubble setup 

The 2D numerical simulations presented in this work have been 
performed with the FLASH code (v. 2.3: lFrvxell et al. 2000), an 
hydrodynamical, adaptive mesh refinement (AMR) code espe- 
cially designed for the study of thermonuclear flashe s. The simu- 
lation setup is very similar to the one described in lapi chino et al.l 
(2006), and we refer the reader to that work for a more detailed 
account. The most relevant differences with respect to it are pre- 
sented in the next sections. A 2D bubble is initialized as a tem- 
perature perturbation, in pressure equilibrium with a background 
medium, in a computational domain with size 5 x 20 km. The 
simulations have Cartesian geometry and were set with a root 
grid of [ 1 6 X 64] grid cells and four additional AMR levels, lead- 
ing to an efi'ective grid size of [256 x 1024], and to an eff'ective 
spatial resolution of 2 x 10"' cm. The initial bubble diameter is 
D = 1 km, considered to be a good order-of-magnitude estimate 
of the size of th e temperature p erturbations in the turbulent WD 
convective flow dWooslev et al.ll2b04.) . The initial distance of the 
bubble from the WD center is set to ^ = 100 km, and the initial 
temperature contrast AT/T is set to 0.07 at the bubble location. 
Both v alues are in the range predicted by Wunsch & Woosle^ 
( |2004|) for the generation of temperature fluctuations in the WD 
core, though these authors refer to the contrast to the core tem- 
perature rather than to the local temperature. Note that, in the 
framework of our simplified model, R and AT/T are not com- 
pletely independent parameters: taking a slightly bigger contrast 
at a given R is like considering a bubble that has already trav- 
elled and grown from a slightly smaller distance to the center. 
Therefore, as long as AT/T is relatively small, its value is not a 
critical issue in the model. 



2.2. WD models 



Unli ke [lap ichino et al.l ( 20061) . we adopt here the WD models 
i lLes afi 



from iLes affre as background for the simulations. 

The models are labelled in Table [T] with the number m in the 
first column, where m = 1,3,4, 6 is related to the initial mass 
of the WD: Minh - (0.6 -1- 0.1 x m) Mp. All models are for a 
cooling age of tg - 0.4 Gyr (see lLesafire et al.ll2006l) . In Table 
[T]the quantity T^ is the central temperature of the model, is 
the central density and g is the gravitational acceleration at a 
distance of 103 km from the WD center The burning zone radius 
of the model rb is defined as the radius within which one-half of 
the luminosity is generated. 
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The details of the modeli ng of the accretion phase are re- 
ported in lLesaffre et alJ (1200 6*). The models are taken at the time 
in the accretion history when the burning time scale is equal to 
the convective turnover time sc ale (this corres ponds to the ig- 
nition criterion a = 1 of Lesaffre et al. I I2OO6I Sect. 2.5). The 
chemical composition of the WD is not set as in llapichino et al.l 
(l2006h (X(12C) = X("'0) = 0.5), but is mapped from the WD 
model. We pick the models from Lesaffre et al. (2006) which use 
values from lUmeda et al.l (Il999h as initial conditions before the 
accretion stage. The C/O ratio in the WD core at ignition ends 
up between 0.74 and 0.46, depending on the model. 

The gravitational acceleration g(r) is derived from the WD 
mass profile as 

gir) = j — (1) 

where G is the gravitational constant and M{r) the mass con- 
tained in the spherical shell with radius r. The acceleration is 
applied to the computational domain as a force pointing down- 
wards along the y-axis, constant in space and time. 

2.3. Nuclear burning and carbon reaction rate 

In the initial WD models the abundances of four chemical 
species are rn a pped, i.e. '^C, '^O, ^^Ne and ^"^Mg. As in 
llapichino et al.l (l2006l) . the hydrostatic nu clear burning is fol- 
lowed with the reaction network iso7 (iTimmes et al.l l2000h 
which evolves seven a-isotopes. As ^^Ne is not available in 
the reaction network, this isotope is substituted by ^°Ne, as 
suggested and verified by Dursi & Timmes (2006). We refer to 
ITimmes et al] (l2000h for the nuclear rates in the reaction net- 
work, with the notable exception of the C-burning, which in 
some calculations is modified according to the following pre- 
scriptions; 

ISpillane et al.l (l2007t) performed a study of the carbon burn- 
ing reactions ^^C('^C, a)^°Ne and ^^C('^C, p)^^Na using y-ray 
spectroscopy with a C target with ultra-low hydrogen contamina- 
tion. The explored energy range spanned from 2. 10 to 4.75 MeV. 
The most relevant result is the detection of a new resonance at 
the very low end of the energy range (Er = 2.14 MeV). 

The effect of this new measurement for the C-burning reac- 
tion rate (the sum of the contributions of the a and proton chan- 
nels) as a function of the temp erature is shown in Fig. [T] With 
respect to the reference value of lCaughlan & Fowleii (Il988h . the 
reaction rate is increased up to 20% in the temperature range be- 
tween 7.0 and 9.8 x 10** K. For comparison, the Gamow energy 
for the C-burning reactions is 1.5+0.3 MeV atT - 5x10^ K, and 
2.1 + 0.4 MeV at r = 8 X 10^ K. The new resonance therefore 
affects the reaction rate mostly in the late hydrostatic burning 
stage, just prior to ignition. The data do not provide informa- 
tion about the burning rate at lower temperatures, relevant for 
the "smouldering phase" of the b urning in mass ive CO WDs 
(5 X 10^ K < r < 7 X 10** K; cf. lNomotol[T982h . An extrapo- 
lation at these temperatures rests upon the new ev aluation of the 
non-re sonant contribution to the reaction rate from' Spillane et all 
(120071) . plus general considerations on the unknown resonance 
structure below 2.0 MeV. Assuming that the observed resonance 
patte rn of the reactions co ntinues in a similar way at lower ener- 
gies, ISpillane et alj (l2007h state that their rates at lower temper- 
ature represent the non-resonant, lower limits to the true values 
at £^ 1.5 MeV. 

On the other hand, it is quite clear that the new rates are 
not relevant for the explosive nucleosynthesis of SNe la, which 
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Fig. 1. Ratio (solid line) between the n ew reaction rate of 
Spi llane et alJ (l2007h for C-burning, and the lCaughlan & Fowleil 
( 198^ rate, in a temperature range relevant for the ignition of 
SNe la. The non-resonant behavior at low temperature has been 
extrapolated from the analytic fit of iSpillane et al. (2007,) . which 
refer to a higher energy range. 



occurs in a different burning regime (nuclear statistical equilib- 
rium) and at r > lO** K, where the corresponding experimental 
uncertainties on the astrophysical factor are much better con- 
strained. 

The relevance of the new rate in the range 7xlO**K<r < 
9.8 X 10^ K motivates therefore the investigation of its effect on 
the ignition process in SNe la, which typically occurs at these 
temperatures. 

Due to the restrictions of the a-network, in iso7 the C- 
burning is implemented in such a way that the reaction rate takes 
into account both the a and the proton channel, but the only re- 
action products are ^"Ne and a. Since the a channel is dominant 
in the considered temperature range, this approximation is sen- 
sible. 

It is worth mentioning that, besides the above cited measure- 
ments, in the last decade several works have been devoted to 
the theoretical prediction of the C-burning reaction rate. Some 
of them (Cussons et al. 2002; Itoh et al. 2003) state that resonant 
screening effects could have a dramatic effect (either suppression 
or enhancement, respectiv ely, for the two cite d pa pers) on the re- 
action rate at low energy. iJiang et al] (l2007h and iGasques et al.l 
(I2007h predict that the fusion cross section has a low-energy hin- 
drance in a range of temperat ures and densit i es whi ch include the 
conditions in the WD core. iGasques et al.l (l2007l) further com- 
pute the SN la ignition curve in this case, finding that, with the 
reduced burning rate, the ignition criteria will be fulfilled at a 
larger temperature. All these theoretical expectations will not be 
addressed explicitly in the present work. 

3. Results 

3.1. Ignition in different WD models 

Four bubble simulations have been performed, one for each WD 
background model listed in Table [1] In these runs, the reaction 
rate is set as in Caughlan & Fowler ( 1988). Figure |2] shows the 
time evolution of the maximum bubble temperature in the four 
simulations. The bubble calculations with models 1 and 3 go to 
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Fig. 2. Maximum bubble temperature as a function of time, in 
four bubble simulations with different WD models, reported in 
the legend. The initial bubble temperature is indicated at f = s, 
the other bubble parameters are set as described in Sect. 12. II 



Fig. 3. Maximum bubble temperature as a function of time in two 
bubble simulations, with the same bubble parameters (described 
in the text) but a different prescription for the reaction rate of 
C-burning, as indicated in the legend. 



thermonuclear runaway, whereas in the cases 4 and 6 the bubble 
is fragmented by the hydrodynamical instabilities and its own 
rising motion before reaching the runaway. The typical evolu- 
ti onary timescales are ab out 25% longer than the ones reported 
in llapichino et al.l (l2006h . because of the lower '^C abundance in 
the present work. 

Two main points explain the differences in the bubble simu- 
lations shown above. First, the initial bubble temperature deter- 
mines the nuclear timescale. A colder temperature perturbation 
has a longer ignition timescale t^, according to 



-7- -22 -4 



(2) 



dWooslev etal.ll2004t llapichino etall l2006h and thus it is dis- 
persed more efficiently by its own buoyant rise. 

Secondly, WD models with larger initial mass are colder and 
denser in the center, and have a larger gravitational acceleration 
at the initial bubble location (cf. Table [TJ. This leads to a larger 
buoyant velocity Vb and thus to smaller dispersion timescales 
Tdisp in models 4 and 6, which can be estimated as 



^disp 



D 

Vb 



(3) 



where D is the bubble diameter 

The ignition condition for a temperature perturbation is 
met when its nuclear timescale is smaller than the dispersion 
timescale. In the case of models 4 and 6, the bubble does not 
go to thermonuclear runaway: for this set of bubble parameters, 
the dispersion is faster than the nuclear burning. Suppose that 
we reduce the bubble distance from the WD center. At the same 
time the ignition timescale decreases (because the temperature 
increases) and the dispersion timescale increases (because g de- 
creases towards the center). Both effects make the ignition more 
likely. 

The main result for this part of our study is therefore a de- 
pendence of the ignition features on the accreting history of the 
WD. WDs with a larger initial mass develop a core which is 
denser and cooler, according to Lesaffre et al. (2006). In these 
models we find that, for a given bubble size, the ignition is fa- 
vored at a smaller distance from the WD center than in models 



with a smaller initial mass. This interpretation is further strength- 
ened by the computation of the burning zone radius r\, of the 
WD models (Table[T] fifth column), defined as the radius within 
which one-half of the luminosity is generated ( Woosley et al.l 
2004). This radius can be considered a good approximation for 
the size of the region where the temperature perturbations arise. 
In denser models, r\, is smaller, with a maximum difference be- 
tween models 1 and 6 of 25%. Temperature perturbations are 
thus likely to be generated closer to the WD center. 

3.2. The role of carbon burning 

The increase of the rate of C-burning in ISpillane et alj (l2007h . 
though not dramatic, has a relevant impact on the evolution of 
temperature perturbations, due to the high power dependence 
on T in Eq. |2] This is shown in Fig. [3] which refers to two 
simulations with the same bubble pa rameters but different refer- 
enc es for the C-buming rate, namelv lCaughlan & Fowlerl (Il988h 
and ISpill ane et al.l (l2007i) . The WD background model is the 
same as in lapichir io et al.l (|2006) and the bubble parameters are 
r = 7.7 X 10*^ K, D = 1 km and = 100 km. Not unexpect- 
edly, the bubble simulation with the new reaction rate reaches 
the thermonuclear runaway, the other does not. From an analyt- 
ical viewpoint, this is well described by the equation linking the 
reaction rate S and t„ ( Wooslev et al. 2004): 





dS\ 




dS 


dT\ 








dT 


dt ] 



CpT 
23$ 



(4) 



where Cp is the specific heat at constant pressure. At T = 
8 X 10^ K, S is 1.2 times larger, thus Tn^new ~ 0.8 Tn_oid, and a 
smaller nuclear timescale favors the bubble ignition before its 
fragmentation. 

This is not a surprising result, and has a limited impact on our 
knowledge of the ignition process. Another interesting question 
pertains to the effect of the new rates not only on a temperature 
perturbation, but on the overall WD structure before ignition. 
In other words, can a l arger rate of C-burni ng modify the WD 
models before ignition (iLesaffre et al.ll2006h and, consequently, 
the features of this process? 
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Table 2. Main features of the two additional WD models com- 
puted to investigate the effects of changing the burning rate 
(quantities are same as in Table[TJ. 



Model 


T, [10** K] 


[10" g cm-^] 


g [10" cm s-^] 


'b [km] 


3.1 


7.68 


2.20 


6.23 


150 


3^ 


7.63 


2.08 


5.90 


154 


3.2 


7.59 


2.00 


5.20 


164 



" As for model 3, the stellar evolution code stalls shortly before the 
actual criterion for ignition is met, hence these values result of an ex- 
trapolation from about 5% lower temperature models. However, we do 
not expect the extrapolation errors to change much either the absolute 
or the relative values of the numbers presented in this table. 



In order to investigate this point, we run two ad d itiona l 
WD models with the code used in iLesaffre et al.l (l2006l) 
(FLASH.THE_TORTOIs4II) in the same conditions as in Model 
3, but we varied the adopted r ate of C-burnin g (Tab le [2]l. Model 
3.1 was run with the r ate of Soilla ne et al.l ([2007) instead of 
the ICaughlan & Fowled ( [1988) rate used for Model 3. The ex- 
trapolation of the rate at low temperature is about half the 
ICaughlan & Fowleil (Il988h rate in the low temperature range be- 
tween 10^ K and 5 X 10^ K (see Fig. [TJ. As seen on Fig. |4] 
the central part of the star spends more time at low temperature 
in that case: as a result the ignition curve is shifted to higher 
densities and the star ignites at a slightly higher d ensity (about 
6% higher, see Table |2]i. The 'Sp illane et al.l (l2007b extrapolated 
burning rate is a lower estimate at low temperature, hence we 
can take the ignition density value of model 3.1 as an upper es- 
timate. The gravitational acceleration follows the same trend as 
the central density. On the contrary, the ignition temperature is 
almost unchanged. This reflects the fact that the reaction rate 
is an extremely steep function of the temperature. The burning 
zone radius is almost unchanged as well. 

Model 3.2 was run with twice the Caughlan & Fowler! d 19881) 
rate, in order to check how big changes in the rate, caused by 
possible discoveries of new low-lying resonances, can affect the 
ignition conditions (for example, the Spillane et al. 200 7 rate 
is, at most, only 20% higher than the Caughlan & FowleJjTgH 
rat e). Since the rate of burnin g is now everywhere higher than 
the lCaughlan & Fowled (Il988l) rate, the ignition curve on Fig.|4] 
is shifted to the left hand side and the star ignites at lower density 
(by about 4%, as seen on Table|2]). As in model 3.1, the ignition 
temperature is almost unchanged. 

To conclude this brief investigation, the main effect of 
the rate is rather surprising at low temperature, where it con- 
trols the density at which the growth of the convective core 
starts, and hence determines the future ignition density. The 
ignition temperature is almost unaffected. The results on ig- 
nition temperatures and burning radii qualitatively agree with 
som e calculations done with t he heuristic convection model 
from lWunsch & WooslevI (|2004 and the reaction rates described 
above, for an isotropic flow and neglecting non-linear effects on 
the burning. 

3.3. Numerical tests on the bubble simulations 

The interpretation of the bubble simulations presented in this 
work cannot be separated from an evaluation of the numeri- 
cal i ssues of the setup. As already discussed in llapichino et al.l 
(12006) . the main features of the bubble physics are the fragmen- 




2.0x10" 2.2x10" 2.4x10" 2.6x 
Central density (g.cm~^) 
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Code web page: | http://astro.ens.fr/?lesaffre| 



Fig. 4. Evolution for the central temperature and density of the 
star in Models 3 (solid). Model 3.1 (dashed) and Model 3.2 
(dotted). Symbols denote the birth of the convective core (trian- 
gles) and ignition points according to various criteria: the burn- 
ing time scale is l/22nd of (crosses) l/8th of (sta rs) or equal to 
(diam onds) the convective turnover timesc ale (see lLesaffre et al.l 
l2006h . Although the position of the crosses and stars corresponds 
to actual output models of the stellar evolution code, the dia- 
monds result from an extrapolation involving models at a lower 
temperature. 



tation due to the buoyant motion, and (if the bubble does not 
reach the thermonuclear runaway) the subsequent dispersion by 
heat diffusion. The latter process is expected on length scales 
(10-100 cm) which are one to two orders of magnitude smaller 
than the effective spatial resolution of the performed simulations. 
Nevertheless, the physics which drives the dispersion, namely 
the interplay between nuclear burning and hydrodynamical in- 
stabilities, is correc tly taken into account in the model, as the 
parameter studies in llapichino et alj ( l2006l) show. Therefore one 
can reasonably expect that the numerical dispersion mimics, on 
larger length scales, the effect of the cooling by heat diffusion at 
unresolved scales. 

One can still argue about other oversimplifications in the 
bubble setup. For example, the choice of a spherical bubble 
shape can be debated. Numerical tests were performed with bub- 
bles in stirred backgrounds, where the bubble is soon deformed 
and the growth of instabilities follows a different pattern with re- 
spect to the standard background model at rest (llapichino et al.l 
i2006i) . The results show that, for well-posed choices of the stir- 
ring field, the turbulence-induced disruption (and thus the bubble 
shape) has not a leading effect in the bubble evolution. 

Closely related to the bubble shape is the geometry choice: 
the 2D Cartesian geometry has been preferred to the more natu- 
ral cylindrical, axisymmetric geometry for numerical reasons, 
after having verified that the results on b ubble dispersion are 
very similar in both cases (Ilapichinoll2005h . 

Another numerical issue concerns the discretization of the 
bubble on the Cartesian grid. This procedure may introduce 
errors, potentially able to seed instabilities that ultimately can 
dominate t he evolution. A resolution study has been therefore 
performed (llapichino et al.l 12006) to verify whether the chosen 
level of resolution allows an adequate modelling of the hydrody- 
namical instabilities and the subsequent fragmentation. In Fig.|5] 
one can notice that different maximum levels of refinement re- 
solve more and more Rayleigh-Taylor unstable modes. On the 



6 



L. lapichino and P. Lesaffre: Uncertainties and robustness of the ignition in SNe la 



































m 


— 












1 




m 


























































































Fig. 5. Comparison of three temperature plots of buoyant bub- 
bles, after a simulation time of 0.5 s, in runs with different max- 
imum level of refinement. The plots enclose only a small part 
of the computational domain (3 x 2.5 km), centered on the bub- 
ble. The choice of the initial parameters is identical in the three 
cases: temperature 7.7 x 10^ K, diameter 1 km, distance from 
the cent er 100 km, and the background WD model is the same 
as in l apichino et al.l ( 120061) . Every mesh element in the plots is 
resolved in 8 x 8 grid cells. The color scale indicates tempera- 
ture, in units of 10^ K. The panels refer to runs with four (upper 
left), five (upper right; this is the reference resolution of the runs 
presented in this work) and six (lower left) AMR levels. 



other hand, the resolution studies of the temperature evolution 
and the bubble dispersion (Figs. 3 and 4 of lapichino et al. 200^ 
clearly have convergent results for the two more resolved runs, 
thus indicating that the bubble evolution is dominated by large- 
scale instabilities, rather than by discretization errors. 



4. Discussion and conclusions 

In this work we explored two previously unaccounted parame- 
ters of the ignition theory of SNe la, and studied how they can 
affect the features of the ignition process. The tools for this study 
are 2D numerical simulations of temperature perturbations in the 
WD core. In th is framework, the bubble model developed by 
llapichino et al. (2006) has several intrinsic and numerical limi- 
tations in its predictive power of the ignition properties of SNe 
la, but it can be profitably used as a "probe" of the influence on 
the ignition process of other physical parameter^ 

Four WD models, derived by the stellar evolution and accre- 
tion history of the progenitors (iLesaffre et al.l l2006). were used 
as background for the bubble evolution. The thermonuclear run- 
away is not reached in the coolest and densest WD models 4 and 
6, resulting from the evolution of initially more massive progeni- 
tors. We infer that, in these models, the ignition region lies closer 
to the WD center, where the gravitational acceleration is smaller 
and the bubble dispersion less effective. 

The simplified bubble model that we used does not allow 
an accurate determination of the extent of the ignition zone. An 



^ For example, following the same approach, lapichin o et al.l ( |2008|) 
propose that the WD composition affects also the ignition features (as 
sug gested bylP ursi & Timmes 2006) and not only the SNe la ex plosion 
(cf. lrimmes eral..,2003. : .Ropke et al.ir2006j : lTownslev et al.ll2009l1 . 



estimate of this quantity is provided by the size of the burning 
zone radius r\,, where one-half of the nuclear energy of the WD 
is produced and where the temperature perturbations are likely 
to be generated. The burning zone radius is in the range between 
170 and 120 km and decreases in more massive progenitors, thus 
supporting the results of bubble simulations. 

Clearly these estimates totally neglect any superimposed 
feature of the WD convection, such as a dipole. In large- 
scale simulations of the progenitor's hydrostatic evolution, on 
the other hand, the spatial resolution is insufficient to resolve 
the generation and ignition (or fragmentation) of the temper- 
ature perturbations. Therefore, we suggest to implement the 
key properties of these objec ts using a stochastic approach 
(cf. lSchmidt & Niemeveil2006l) . 

The second question addressed by this work is about the nu- 
clear physics uncertainti es in the ignition th eory. Starting from 
the new measurement of lSpillane et al.l (|2007|) . we study the im- 
pact of a different (larger) reaction rate for the C-burning on the 
ignition features in SNe la. This analysis is not limited to the 
temperature perturbations, but is extended to the WD structure 
by computing new stellar evolution tracks with various rates of 
C-burning. According to these models, the ignition features only 
have a mild dependence on the rate at low temperature (where 
lower rates favor higher ignition densities) and almost no de- 
pendence on the rate at high temperature. A somewhat similar 
analysis has been performed by Cooper et al. (2009): they study 
the possible effect of a hypothetical resonance at £ = 1 .5 MeV 
by constructing ignition curves for the progenitor They also 
find that such a feature would play no major role on the igni- 
tion process, though lower-lying resonances could. It is impor- 
tant to stress that the effect of different ignition conditions is 
mani fold and not merely related to the location of the ignition 
zone: ICooper et al] (12009 ^. for example, c onjectu re an impact on 
the iron-peak explosiv e nucleosynthesis dPiro & Bildsten 200^ 
Chamulak et al. 2008^ and on the velocity of the flame front 
(Ropke et al. 2006a ). 

The whole field of the theoretical modeling of SNe la is 
reaching in recent years a mature stage, with growing interest- 
ing insights in the microphys ics of these events (like fo r exam- 
ple the mechanism for DDT, Ropke & Niemever 2007; Ropk^ 
12007; Wooslev..2007.) . Moreover, modern simulations provide 
a more and more deta i led comparisons with the observables 
(•Hi llebrandt et al.ll2007t iMazzaU et al.ll2007t iKasen & WooslevI 
200^" 

The use of SNe la for precision cosmology has fi- 
nally come on safer grounds, thanks to a deeper the oretical un- 
derstanding of the involv ed physical parameters (iRopke et al.l 
l2007d: iKasen etal]|2009h . At the same time, novel low-Mach 
number schemes permit to fill the gap between the evolution 
of the proge nitor and the onset of the explosion (Zingale et al] 
I200d l2()09l) . In this context, the present study on the ignition 
parameters, although performed with very simplified tools, can 
drive further and more sophisticated analysis towards interesting 
issues, which call for a deeper investigation. 
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